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2 ■ ABSTRACT 

' We calculate the major pair fraction and derive the major merger fraction and rate for 82 

massive (M* > 10 n M©) galaxies at 1.7 < z < 3.0 utilising deep HST NICMOS data taken 
in the GOODS North and South fields. For the first time, our NICMOS data provides imaging 
with sufficient angular resolution and depth to collate a sufficiently large sample of massive 
galaxies at z > 1 .5 to reliably measure their pair fraction history. We find strong evidence that 
\^ ■ the pair fraction of massive galaxies evolves with redshift. We calculate a pair fraction of f m = 

| 0.29 +/- 0.06 for our whole sample at 1.7 < z < 3.0. Specifically, we fit a power law function 

Q\ . of the form f m — /o(l + z) m to a combined sample of low redshift data from Conselice et 

al. (2007) and recently acquired high redshift data from the GOODS NICMOS Survey. We 
find a best fit to the free parameters of f = 0.008 +/- 0.003 and m = 3.0 +/- 0.4. We go on 
to fit a theoretically motivated Press-Schechter curve to this data. This Press-Schechter fit, 
qq ' and the data, show no sign of levelling off or turning over, implying that the merger fraction 

(~) of massive galaxies continues to rise with redshift out to z ~ 3. Since previous work has 

established that the merger fraction for lower mass galaxies turns over at z ~ 1.5 - 2.0, this 
is evidence that higher mass galaxies experience more mergers earlier than their lower mass 
counterparts, i.e. a galaxy assembly downsizing. Finally, we calculate a merger rate at z = 2.6 
5_i ■ of < 5 x 10 5 Gpc~ 3 Gyr -1 , which experiences no significant change to < 1.2 x 10 5 

Gpc~ 3 Gyr -1 at z = 0.5. This corresponds to an average M* > 10 u Mq galaxy experiencing 
1.7 +/- 0.5 mergers between z = 3 and z = 0. 

Key words: galaxies: formation, galaxies: mergers, galaxies: evolution, galaxies: massive, 
cosmology 
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1 INTRODUCTION 

Hierarchical assembly is the established leading theory for the evo- 
lution of structure in the universe from the viewpoint of Cold Dark 
Matter (CDM) theory. This states that larger structures form from 
the merging of smaller structures. Thus CDM halos grow in size, 
with the largest structures forming latest in the history of the uni- 
verse. Galaxies are often thought to form principally in line with 
halo mergers. If they do, it is likely that merging of smaller galaxies 
to create more massive ones is the major factor in the evolution of 
galaxies over cosmic time. Alternatives to a merger history of mas- 
sive galaxies include rapid collapse mechanisms, whereby galaxies 
form over very short time-scales in the early universe, and then do 
not significantly merge during the rest of their lifetime. 



The aim of this paper is to address the question: what mech- 
anism drives galaxy evolution with regards to the most massive 
(M» > 10 Mq) galaxies in the universe? We can observation- 
ally test this by calculating the merger fraction of massive galaxies 
at different redshifts. Previous work by Conselice et al. (2007) ex- 
plores this problem for similar mass galaxies out to z ~ 1 .4, using 
morphological techniques. We extend this work using a close-pair 
method to z ~ 3, using data from the GOODS NICMOS Survey. 

There are several different ways to locate merging galaxies 
which roughly relate to the stage of the merger. The most direct 
method is to look at morphological disturbances in galaxies, e.g. 
Conselice et al. (2003, 2006 and 2008). In this approach one selects 
galaxies with asymmetries, or distortions in their morphologies, 
above a certain threshold, and define these as merging systems. 
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This can be carried out by eye, or via computational methods, e.g. 
CAS, Gini and M20 parameters (see Conselice 2003, 2006, Con- 
selice at al. 2008 and Lotz et al. 2008a). Alternatively, one could 
also look in principle for secondary features such as an abrupt in- 
crease in star formation rate. This is because the merging of gas 
rich galaxies will often give rise to an increase in star formation, 
and hence a distinctive spectral signature, see e.g. Lin et al. (2008) 
and Mathis et al. (2005). 

Also one can look for galaxies in close proximity and count 
the number of galaxies in apparent pairs in a given sample. This 
pair counting method requires much less angular resolution than 
the morphological approaches, and can be extended to higher red- 
shift more directly than the others. The close pair method does not 
precisely trace merging, however, but rather looks for likely po- 
tential future mergers. Thus, to go from a close pair fraction to a 
merger fraction requires that assumptions be made. Ultimately, the 
pair count method measures mergers before they happen, the mor- 
phological approach measures merging while it happens, and meth- 
ods involving increased star formation rates trace merging during 
and after the event. 

Merging is already known to be of paramount importance in 
driving the rate of star formation, the formation and evolution of su- 
permassive black holes, and as a mechanism for increasing galaxy 
mass. Further, mergers can be used to trace the formation history of 
galaxies, probing where, when, and how galaxies form. Conselice 
et al. (2007) and Rawat et al. (2008) find that massive galaxies 
(M, > 10 n M Q ) out to z ~ 1.5 have an increasing merger fraction 
with redshift. It has been noted for lower mass systems that there 
is a levelling off and eventual turn around in the merger fraction at 
z < 2 (see Conselice et al. 2008). In order to probe whether the 
same levelling off and turn around in the merger fraction occurs 
for M* > 10 11 Mq galaxies, it is necessary to explore the merger 
fraction of these massive galaxies at higher redshift. 

Consequently, we have collated a sample of 82 galaxies with 
M, > 10 n M Q , selected from the GOODS North and GOODS 
South fields. These are imaged as part of the GOODS NICMOS 
Survey (GNS). The GOODS field has the multiband coverage 
needed to estimate good photomentric redshifts, in addition the 
depth and angular resolution from NICMOS allows us to resolve 
significantly fainter objects around the primary targets to distances 
of a few kpc. Thus, for the first time, we can probe the pair frac- 
tion history to z ~ 3. We argue that the increase in merger fraction 
with redshift continues for massive galaxies in our sample out to 
z ~ 3. The merger fraction shows no sign of levelling off or turn- 
ing around at 2 < 2 < 3, unlike for lower mass systems observed 
by Conselice et al. (2008). Thus, we find that, complimentary to 
the star formation results of Bundy et al. (2006) and Pipino et al. 
(2008), the most massive galaxies appear to form by mergers at a 
higher rate in the early universe than less massive ones. 

Throughout this paper we assume a ACDM Cosmology with: 
Ho = 70 km s _1 Mpc -1 , VL m = 0.3, S1a = 0.7, and adopt AB mag- 
nitude units. 



2 DATA AND OBSERVATIONS 

The GOODS NICMOS Survey imaged a total of ~ 8000 galaxies 
in the F160W (H) band, utilising 180 orbits and 60 pointings of the 
HST NICMOS-3 camera. These pointings are centred around mas- 
sive galaxies at z = 1.7 - 3 at 3 orbits depth. Each tile (51.2"x51.2", 
0.203"/pix) was observed in 6 exposures that combine to form im- 
ages with a pixel scale of 0.1", and a point spread function (PSF) of 



~ 0.3" full width half maximum (FWHM). See Magee, Bowens & 
Illingworth (2007) for details relating to the data reduction proce- 
dure. The pointings were chosen and optimised to contain as many 
high-mass (M, > 10 11 M Q ) galaxies as possible. The selection of 
these is outlined in Conselice et al. (2008, in prep.). These galax- 
ies consist of high redshift galaxies selected by various optical- 
to-infrared colour techniques (see Papovich et al. 2006, Yan et al. 
2004 and Daddi et al. 2007). A total of 82 galaxies were found with 
M» > 1O 11 M , with photometric and spectroscopic redshifts be- 
tween z = 1.7 and z = 3. Limiting magnitudes reached are H = 26.5 
(10a). 

There is a wealth of observational data covering the GOODS 
field. As such, the masses and photometric redshifts of our sample 
of massive galaxies are calculated using data from the U-band to 
the H-band (e.g. Giavalisco et al. 2004). The stellar masses were 
measured utilising standard multi-colour stellar population fitting 
techniques, producing uncertainties of ~ 0.2 dex. The details of this 
procedure can be found in, e.g. Bundy et al. (2006) and Conselice et 
al. (2007). The stellar masses were calculated assuming a Chabrier 
initial mass function and by producing model spectral energy dis- 
tributions (SEDs) constructed by Bruzel & Chariot (2003) via stel- 
lar population synthesis models parametrised by an exponentially 
declining star formation history. The model SEDs are then fit to 
the observed SEDs of each galaxy to obtain a stellar mass. Recent 
work by Conroy, Gunn and White (2008), however, has suggested 
that current methods for measuring stellar masses may have greater 
errors than previously thought (~ 0.3 dex). Although newer meth- 
ods have been developed utilising different stellar evolution tracks, 
we find that these newer models do not significantly affect our mea- 
sured stellar masses, and our stellar masses are accurate to within 
a factor of a few. For a discussion of this see Trujillo et al. (2007) 
and Conselice et al. (2008). Due to the steepness of the stellar-mass 
function at high masses, Poisson errors in measuring the mass of 
galaxies may lead to more lower mass galaxies being counted in 
our high mass sample. To model this we ran a Monte-Carlo sim- 
ulation based on the stellar mass function for high z galaxies in 
Fontana et al. (2006). We found that ~ 8% of lower mass galaxies 
systematically infiltrated our higher mass sample. This effect is rel- 
atively small and will not significantly affect our results. Moreover, 
since we go on to note a dramatic increase in the pair fraction of 
massive galaxies, when compared to less massive ones, this effect 
cannot be responsible because contamination of lower mass sys- 
tems would result in the pair fraction of higher mass galaxies being 
lower (see e.g. Conselice et al. 2007). 

Our photometric redshifts are determined via standard tech- 
niques (e.g. Conselice et al. 2008). These contain additional uncer- 
tainty on top of the stellar mass issues. Additionally, we find seven 
spectroscopic redshifts from the literature for our sample. Using 
the GOODS/VIMOS DR1 (see Popesso et al. 2008) we find three 
matches with = 0.026, and four spectroscopic redshifts from 
a compilation of redshifts from the literature (see Wuyts et al. 2008) 
giving fa = 0.034. 

We use SExtractor to make catalogs for all detected galaxies in 
the GNS, removing stars and spurious detections. We then carefully 
check by eye that the deblending is accurate and that we are not ar- 
tificially merging separate galaxies or splitting individual galaxies 
into parts. For the massive galaxies studied in this paper, we put in 
particular effort to ensure that the sources are extracted accurately. 
Further, the choice of studying major mergers with mass ratios of 
at least 1:4 of the primary source (+/- 1.5 mag.) allowed us to be 
complete. 
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3 RESULTS 

3.1 Merger Fraction 

3.1.1 Close Pair Method 

We follow an approach similar to Patton et al. (2000) for calculating 
major merger fractions for our galaxies based on the pair method. 
First we separate our galaxies into two redshift bins, with 1.7 < 
z < 2.3 for the first group (containing 44 galaxies) and 2.3 < 
z < 3.0 for the second group (containing 38 galaxies). The method 
for assigning potential pairs is relatively straightforward. We assign 
any galaxy within 30 kpc, in physical units (assuming Ho = 70 km 
s _1 Mpc -1 ), of our host galaxy to be a count, if it is within +/- 
1.5 of the host galaxy magnitude. To calculate a pair fraction we 
sum up the number of galaxies within 30 kpc of all of our host 
galaxies and divide by the total number of galaxies. This would be 
an accurate pair fraction if we had precise redshifts for all of the 
galaxies in our sample, and included in the summation only those 
at the same redshift as the host. As we have redshifts only for our 
host galaxies, we do not know if any particular galaxy within 30 
kpc is a real pair or just foreground or background contamination. 

We correct for contamination by calculating the probability 
of a close galaxy being a pair by chance from our surface number 
counts. These number counts are taken around the objects to min- 
imise clustering effects. We then subtract this correction from our 
running pair total. Specifically we calculate: 

/.m+1.5 

corr = / p(m') x 7r(r| ofepc - rl kpc )dm' (1) 

Jm-1.5 

Where p(m') is the surface density of galaxies in the survey, in the 
magnitude range +/- 1.5 of the magnitude of the host galaxy, r^o^pc 
and r 5 k P c are distances of 30 and 5 kpc from the host galaxy's cen- 
tre, respectively. Essentially this expression corresponds to calcu- 
lating how many galaxies in the whole survey are within +/- 1.5 of 
the magnitude of the host galaxy, and dividing this number by the 
total area of the survey. We then multiply by the area within the an- 
nulus marked out by the radii 5 and 30 kpc, which is done to avoid 
miscounting due to blending with the host galaxy. This gives an ex- 
pectation value for the number of galaxies one would expect to see 
within 30 kpc of the host galaxy purely by chance. Although mas- 
sive galaxies at high redshifts are highly clustered, our correction is 
based upon the densities of objects centred around these systems. 
This value is then subtracted from the number counts to obtain a 
major pair fraction thus: 

- i=N 

f m ~Tf ( counts i - corr (2) 

Where N is the total number of galaxies in the summation. We use 
the values of +/- 1.5 for the magnitude range to select major (1:4) 
mergers only (where we are complete). Furthermore, we adopt the 
convention of setting 30 kpc as our pair distance to be in line with 
Patton et al. (2000) and Bundy et al. (2004). This allows us to make 
fruitful comparisons, and also follows rough theoretical arguments 
for the likelihood of a major close pair becoming a major merger 
in a short (~ 400 Myr) time-scale. The merger fractions derived 
via close pair methods at d = 30 kpc also have a similar timescale 
to the merger fractions calculated via CAS methods. For further 
discussion on this and the characteristic timescales relevant to these 
two approaches, see Conselice, Yang & Bluck (2008). 

We find a pair fraction of f m = 0.29 +/- 0.06 for the whole 
sample at 1.7 < z < 3.0. Within the range 1.7 < z < 2.3 we 
calculate a pair fraction f m =0.19 +/- 0.07, and within the range 
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Figure 1. The merger fraction evolution of M» > 10 11 Mq galaxies. The 
red sqaures are taken from POWIR data, with merger fractions calculated 
via CAS morphologies (Conselice et al. 2007). The blue circles are taken 
from POWIR data, with merger fractions calulated via close-pair methods. 
Blue triangles are taken from the GNS, with merger fractions calculated 
via close-pair methods. The green circle represents the local universe value 
(calculated in de Propris et al. 2007). Solid line is a best fit power law, to 
the high z data, of the form f m = 0.008(1 + z) , with dotted line being 
a best fit power law exponential, to the high z data, of the form / m = 
0.008(1 + z) ' 3 exp(1.0(l + z) 2 ). 

2.3 < z < 3.0 we calculate a pair fraction f m = 0.40 +/- 0.10. 
These are considerably higher values than those found by Conselice 
et al. (2008) for similar mass objects at z < 1.4. We do not make 
any correction to transform our pair fractions into merger fractions, 
unlike Patton et al. (2000) and Rawat et al. (2008). This is because 
both morphological and close pair methods likely trace merging but 
have different timescales. A detailed investigation on the differing 
timescales traced by morphological and pair methods in Conselice, 
Yang & Bluck (2008) suggests that the timescale for 30 kpc pairs 
and morphologically selected mergers is very close. Moreover, we 
go on to calaculate the pair fractions for the POWIR survey mas- 
sive galaxies, to directly compare these values to the GNS massive 
galaxies (see Fig. 1). 

N-body simulations from Wetzel et al. (2008) suggest that 
pair fraction methods may underestimate the number of true ma- 
jor mergers, as pairs at higher separations may also merge. If this 
is true, it suggests that the real merger fractions may be higher 
than what we calculate. But, since the calculated merger fractions 
in this paper are very high, this will not change the thrust of our 
conclusions and may even make them stronger. Our pair fractions 
are plotted alongside lower redshift points from the POWIR Sur- 
vey, with merger fractions estimated via CAS methods (Conselice 
et al. 2007) and via close pairs (calculated in this paper), and with 
a local universe value calculated in de Propis et al. (2007) based 
on morphological methods in Fig. 1. For a detailed explanation of 
morphological techniques see Conselice et al. (2007). 

3.1.2 Merger Fraction Evolution 

We plot a fit of the form f m = / (1 + z) m to all POWIR and 
GNS points (solid line in Fig. 1), and find a best fit to the free 
parameters of: f = 0.008 +/- 0.003 and m = 3.0 +/- 0.4. Our value 
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of fo compares very favourably to the accepted local universe value 
of the merger fraction for massive galaxies of 0.009 +/- 0.002 at z 
= 0.05 (de Propris et al. 2007). We also fit a power law exponential 
curve of the form/ m = a(l+z) b exp(c(l+^) 2 ) to all POWIR and 
GNS points (dotted line in Fig. 1) (Carlberg, 1990). We find a best 
fit to the free parameters of: a = 0.008 +/- 0.002, b = 0.3 +/- 0.2 and 
c = 1.0 +/- 0.6. There is very poor agreement, however, between the 
power law exponential prediction for the local universe value and 
the de Propris et al. (2007) point. This parameterisation also shows 
no sign of turning over at high z, unlike for the lower mass systems 
in Conselice et al. (2008). As such, this indicates that the merger 
fraction continues to increase with redshift out to z ~ 3, implying 
that higher mass objects are experiencing relatively more mergers 
than their lower mass counterparts at high redshift, compared to the 
lower mass data in Conselice et al. (2008). Further, it is also evident 
that the simple power law curve fits our data better than the more 
complex power law exponential. 
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3.2 Merger Rates 

In order to calculate merger rates from our merger fractions we 
must assume a time-scale (i~ m ) over which merging is occuring for 
galaxies in a major pair. Adopting the simulation results found in 
Lotz et al. (2008b), we take r m = 0.4 +/- 0.2 Gyr for close pairs 
at d = 30 kpc, and r m = 1.0 +/- 0.2 Gyr for morphological CAS 
measurements. We also need to know the comoving number den- 
sity, n(z), for Mi. > 10 11 Mq galaxies at different redshifts. We 
calculate this using data from Drory (2005) for the GOODS South 
field. The typical errors of these values are +/- 20 % of the value. 
It is also important to use the galaxy merger fraction (f gm ) as op- 
posed to the merger fraction calculated above. The galaxy merger 
fraction relates to the merger fraction (f m ) by: 

2 X fm 



fgm — 



(3) 



l+fm 

This gives the number of galaxies merging as opposed to the num- 
ber of mergers (/ m ), divided by the number of hosts in the sample. 
For a more detailed derivation of this and discussion on its applica- 
tion see Conselice (2006). We calculate the merger rate thus: 



K 0) = fgm(z)n(z)r n 



(4) 



where n(z) is the comoving number density of galaxies and 
fgm(z) is the galaxy merger fraction (calculated above). We calcu- 
late this rate at z = 2.6 as K < 5 x 10 5 Gpc" 3 Gyr" 1 and at z = 0.5 
as 5ft < 1.2 x 10 s Gpc -3 Gyr -1 . These rate calculations are, how- 
ever, notoriously difficult to perform accurately as there are signifi- 
cant errors associated with the calculation of the comoving number 
density (n), the galaxy merger fraction (fgm), and the characteristic 
timescale for a merging system being detected (r m ). When added in 
quadrature, the resultant final error on the merger rate is very large. 
Therefore, unlike for the merger fraction, the major merger rate is 
consistent with having no redshift evolution. There are two com- 
peting variables here: the comoving number density, n(z), which 
decreases with increasing redshift, and the galaxy merger fraction, 
fgm(z), which increases with increasing redshift. These appear to 
combine to roughly negate any redshift evolution of the merger rate 
R(z) (Fig. 2). 

A useful term may be defined (r) which is essentially the char- 
acteristic time between mergers that an average galaxy experiences 
at a given redshift. It is defind as: 



r 



Tn 

fgr 



(5) 



Figure 2. The major merger rate, !R(z). Red squares are taken from POWIR 
data with merger fractions calculated via morphological CAS methods, blue 
circles are taken from POWIR data with merger fractions calculated via 
close pair methods and the blue triangles are high redshift data taken from 
the GNS with merger fractions derived from close pair methods. 
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Figure 3. The redshift evolution of the average time between mergers 
r = T m / fgm - Red squares are taken from POWIR data with merger frac- 
tions calculated via morphological CAS methods, blue circles are taken 
from POWIR data with merger fractions calculated via close pair meth- 
ods and the blue triangles are high redshift data taken from the GNS with 
merger fractions derived from close pair methods. The solid line is the fit 

r = 12(1 + z)- 1 - 6 . 



This is plotted in Fig. 3. It is clear that F evolves with redshift from 
~ 1.5 Gyr at z = 3 to ~ 12 Gyr at z = 0. Further, we fit a curve of 
the form T = To(l + z)~ p , with a best fit to the free parameters 
of: To = 12 +/- 3 Gyr and p = 1.6 + /- 0.4. This represents a steep 
increase in the time taken for galaxies to merge. We calculate the 
number of expected mergers (N m ) for a given galaxy between two 
redshifts using this power law fit. Specifically we calculate: 



N m = 



1 



■dt = 



tn 



dz 



T(z) (1 + z) E(z) 



(6) 
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where T(z) is the characteristic time between mergers, tn is the 
Hubble time, and the parameter E(z) — [Q,m(1 + z) 3 + f2fc(l + 
zf + Q A r 1/2 = H' 1 (z). Calculating this from z = 3 to z = 0, 
we obtained a value of N m = 1.7 +/- 0.5 major mergers per galaxy, 
with T m = 0.4 Gyr. 



4 DISCUSSION 

Conselice et al. (2008) find that galaxies with masses in the range 
8 < log(M«/M ) < 9 and 9 < log(M*/Af Q ) < 10 have apeak 
in their merger fraction at around z = 1.5, with galaxies with stel- 
lar masses log(Af„/Af0) > 10 appearing to peak in their merger 
fraction later, at around z = 2. In this paper we have found that for 
galaxies with log(M„/M0) > 11 the merger fraction continues to 
increase out to z ~ 3, indicating that the peak must occur at z > 
3. This provides further evidence that more massive galaxies un- 
dergo a greater number of major mergers earlier than less massive 
systems. 

This method of calculating pair fractions is based on not- 
ing overdensities of galaxies around host galaxies, based on back- 
ground counts. However, the clustering of galaxies increases with 
redshift (see Bell et al. 2006). So, this effect may be partly ex- 
plained by the correlation function of galaxies increasing with red- 
shift. Since the parameters of this function are poorly known obser- 
vationally, no attempt has been made to correct for this issue. Con- 
sequently, our pair fraction results from both the GNS and POWIR 
surveys may be interpreted as the 2 point correlation at 30 kpc for 
these massive galaxies. This clearly evolves and it is likely true that 
the actual pair fraction evolves in line with this correlation. Con- 
tamination from line of sight projection at scales larger than 30 
kpc will be minimised by the fact that we calculate our correction 
around the galaxies being measured. Further, we have recently ac- 
quired morphological CAS measurements on the GNS galaxies as 
well as the POWIR ones. This gives a consistent merger fraction 
(/m ~ 0.3) to the pair methods. We explore this in a forthcoming 
paper, and elude to the relation between morphologiocal and pair 
methods in Conselice, Yang & Bluck (2008). 

Despite the large errors associated with calculating merger 
rates, we have demonstrated that there is little or no evolution in the 
rate of major mergers with redshift. The maximum evolution per- 
mitted by the errors would only amount to increasing the merger 
rate by a factor of a few. Of greater significance, however, is our 
calculation of the redshift dependence on the characteristic time 
between mergers, T(z). This quantity decreases significantly with 
redshift. As such, in the early universe (z ~ 3) there was a much 
shorter time between mergers than there is today. 



5 SUMMARY 

We derive the major merger fractions for massive (M, > 10 n M Q ) 
galaxies in the GNS at two redshifts. We find a merger fraction 
of f m = 0.19 +/- 0.07 at z = 2.0 and f m = 0.40 +/- 0.10 at z = 
2.6. This indicates, when compared to data from Conselice et al. 
(2007), De Propis et al. (2007), and POWIR pair fractions (also 
calculated in this paper) in Fig. 2, that there is a strong correlation 
between merger fraction and redshift for the massive galaxies in our 
sample. Furthermore, we fit a function of the form f m = /o(l + 
z) m to our data, and the data from Conselice et al. (2007) at lower 
redshifts for M* > 10 n M Q galaxies, finding a best fit to the two 
free parameters of: / = 0.008 +/- 0.003 and m = 3.0 +/- 0.4. This 



indicates that the merger fraction for massive galaxies continues to 
increase out to z ~ 3. This implies that more massive galaxies have 
higher numbers of mergers than less massive ones, since pevious 
work has found that the merger fraction for lower mass galaxies 
levels off and declines with increasing redshift in the range z ~ 1.5 
- 2.0 (Conselice et al. 2008). 

We calculate the major merger rate finding no significant evo- 
lution from K < 5 x 10 5 Gpc" 3 Gyr" 1 at z = 2.6 to ft < 1.2 x 10 5 
Gpc -3 Gyr -1 at z = 0.5. Moreover, a steep evolution is ruled out. 
We also calculate the evolution of the characteristic time between 
mergers T(z), finding a rapid decrease from ~ 12 Gyr at z = to 
~ 1.5 Gyr at z = 3. This indicates that massive galaxies take longer 
between individual mergers in the local universe than at high z, but 
that there are more massive galaxies now than there were in the 
early universe, so, the rate remains constant. By integrating over 
the best fit to this curve, we calculate the average number of ma- 
jor mergers a massive (M» > 10 11 Mq) galaxy would experience 
between z = 3 and z = is N m = 1.7 +/- 0.5. 

We would like to thank our collaborators on the GNS team, 
particularly Fernando Bruitrago and Ignacio Trujillo, for their work 
on this project and help in producing this paper, and Samantha 
Penny for her assistance. We also gratefully acknowledge support 
from the STFC. 
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